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ABSTRACT 

Formalisms for both non-relativistic as well as relativistic versions of field emission of 
electrons in presence of strong quantizing magnetic field, relevant for strongly magne- 
| tized neutron stars or magnetars are developed. In the non-relativistic scenario, where 

■ electrons obey Schrodinger equation, we have noticed that when Landau levels are 

populated for electrons in presence of strong quantizing magnetic field the transmis- 
sion probability exactly vanishes unless the electrons are spin polarized in the opposite 
direction to the external magnetic field. On the other hand, the cold electron emission 
under the influence of strong electrostatic field at the poles is totally forbidden from 
t— > ' the surface of those compact objects for which the surface magnetic field strength is 

3> I0 15 G (in the eventuality that they may exist). Whereas in the relativistic case, 
where the electrons obey Dirac equation, the presence of strong quantizing magnetic 
field completely forbids the emission of electrons from the surface of compact objects 
if B > 10 13 G. 
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O ' i INTRODUCTION 

There are mainly three kinds of electron emission processes from metal surface, they are: (i) thermionic emission, (ii) photo- 
' electric emission and (iii) cold emission or field emission. 

The field emission or cold emission, which we have investigated in the present article in the context of strongly magnetized 
rS ' neutron stars or magnetars, is an electron emission process induced by strong external electrostatic field at zero or extremely 
, low temperature. Field emission can happen from solid and liquid surfaces, or from individual atoms. It has been noticed that 
the field emission from metals occurs in presence of high electric field: the gradients are typically higher than 1000 volts per 
micron and the emission is strongly dependent upon the work function of the material. Unlike the thermionic emission and 
photo-emission of electrons, the field emission process can only be explained by quantum tunneling of electrons, which has 
no counter classical explanation. However, for general type surface barrier, this purely quantum mechanical problem can not 
be solved exactly, a semi- classical approximation, known as WKB (the name is an acronym for Wentzel-Kramers-Brilloun) 
is needed to get tunneling coefficients. Now to explain cold electron emission from metals, one may assume that because 
of quantum fluctuation, electrons from the sea of conduction electrons (degenerate electron gas) always try to tunnel out 
through the metallic surface (surface barrier). However, as soon as an electron comes out, it induces an image charge on the 
metal surface, which pulls it back and does not allow this emitted electron to move far away from metal surface in the atomic 
scale. But if some strong attractive electrostatic field is applied near the metallic surface, then depending on the Fermi energy 
of electrons, the height of the surface barrier and the local work function, the electrons may overcome the effect of image 
charge and get liberated. Since the external strong electric field is causing such emission and does not depend on the thermal 
properties of the metal, even the metal can be at zero temperature, it is called field emission or cold emission. 

The theory of field emission from bulk metals was first pro posed by Fowler and Nordh eim in an epoch makin g pape r 



in the proceedings of Royal Society of London in the year 1928 iFowler fc Nordheim I (|l928h (see also IStern et al. 1 (|l929h 
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Jensen I (| 19951 ); iForbes fc Deanel (|2007t ) ; iLiang fc Chen I (|2008l ) for further discussion). Fowler-Nordheim tunneling is the 



wave-mechanical tunneling of electrons through a triangular type barrier produced at the surface of an electron conductor by 
applying a very high electric field. 

Now the cold emission or field emission processes not only have significance in the terrestrial laboratories, but is found to 
be equally important in the electron emission processes from cold and co mpact stellar objects, such as neutron stars. In the 
case of a rotating neutron star, the existing large mag net ic field, > 10 12 G 
radio pulsars or > 10 15 G for the magnetars inner field iMereghetti 
approximately given by F ~ 2 x 10 8 P -1 -Bi2 Volt cm -1 and is parallel to 



Shapiro fc Teukolskv I (|l983l ) for the conventional 
produces a strong electric field at the poles, 



B at the poles Shapiro fc Teukolskv ( 1983}), here 
P is the time period of the neutron star in second and B12 is the measure of magnetic field strength in the units of 10 12 . 

At the proximity of polar region of a strongly magnetized neutron star, the potential difference changes almost linearly 
with distance from the polar cap, which is a region very close to the magnetic pole. The repulsive surface barrier in combination 
with this attractive potential, forms a triangular type barrier at the poles. Therefore if electron emission from the poles of 
neutron stars is field emission type, then Fowler-Nordheim equation with proper modification may be used to investigate 
such emission process. Now, the study of plasma formation in a pulsar magneto-sphere is a quite old but still an unresolved 
astro p hysical issue, in partic u lar th e magneto-spheres of strongly magnetized neutron stars /magnetars iMolofeev et al. I (J2004, 



2005); llstomin fc Sobvanin 



(120071 ). In the formation of magneto-spheric plasma, it is generally assumed that there must 
be an initial high energy electron flux from the magnetized neutron stars. Since the magnetic field at the poles of neutron 
stars/magnetars are strong enough, the emitted electrons flow only along the magnetic field lines. The flow of high energy 
electrons along the direction of magnetic lines of forces and their penetration through the light cylinder is conventionally 
pictured with the current carrying conductors. Naturally, if the conductor is broken near the pulsar surface the entire potential 
difference will be developed across a thin gap, called polar gap. This is of course based on the assumption that above a critical 
height from the polar gap, because of high electrical conductivity of the plasma, the electric field F, parallel to the magnetic 
field near the poles is quenched. Further, a steady acceleration of electrons originating at the polar region of neutron stars, 
travelling along the field lines, will produce magnetically convertible curvature 7-rays. If these curvature 7-ray photons have 
energies > 2m e c 2 (with m e is the electron rest mass and c is the velocity of light), then pairs of e~ — e + will be produced 
in enormous amount wit h very high efficiency near th e polar gap. Th e se produced e~ — e + pairs form what is known as the 



magne t o-spheri c plasm a Shapiro fc Teukolskv dl983[); Jessner et al. 1 2001 ); Ruderman fc Sutherland 1 19751) ; Diver et al 



[ 2009h : lMicheT1 1 1981 [20041 ): Irlarding fc Lai I (200rj ): lRuderman I l| 19711 1. The cold emission, therefore plays a significant role 
in magneto-spheric plasma formation. In turn, the motion of charged particles in the magnetosphere in presence of strong 
magnetic field causes pulsar emission in the form of synchrotron radiation. Therefore the cold emission process indirectly also 
affects the intensity of synchrotron radiation 



Further the e xactly solvable models with simple type tunneling barrier lead to equations Fowler fc Nordheim (1928) 



Stern et al. 1(1 19291) that underestimates the emission current density by a factor of 1000 or more. If a more realistic type barrier 



model is used by inserting an exact surface potential in the simplest form of the Schrodinger equation, then a complicated 
mathematical problem arises over the resulting differential equation. It is in principle therefore mathematically impossible to 
solve the equation exactly in terms of the usual functions of mathematical physics, or in any simple way. 

Moreover, to the best of our knowledge, neither the non-relativistic nor the relativistic version of cold emission processes 
in presence of strong quantizing magnetic field, relevant for electron emission from the poles of strongly magnetized neutron 
stars/magnetars, even with simple type potential barriers have been properly investigated. In the conventional pulsar model 
it is generally assumed that the emission of electrons and thereby formation of magnetosphere is mainly caused by strong 
electric field at the polar region which is produced by the strong magnetic field of rotating neutron stars. Taking this physical 
picture into consideration, in this article, we have developed formalisms for both non-relativistic and relativistic scenarios of 
field emissions for electrons from the poles of neutron stars with 10 10 G < B < 10 17 G. 

In the next section we have studied the effect of strong quantizing magnetic field on the field emission of electrons for 
the non-relativistic case. In section-3, we have repeated the same calculation for the relativistic scenario. In the last section 
we have given the conclusions and future prospects of this work. 



2 EFFECT OF STRONG QUANTIZING MAGNETIC FIELD ON COLD EMISSION: 
NON-RELATIVISTIC SCENARIO 

To develop the modified version of field emission process for the electrons in the non-relativistic scenario in presence of 
a s trong quantizing ma g netic field, we have followed the basic calculation presented in the seminal Royal Society paper 
by Fowler fc Nordheim ( 19281 ). In the modified version, we assume a cylindrical co-ordinate system (p,9,z) and the con- 



stant magnetic field B is a long positive z-direction, with the usual gauge for the vector potential A = (B x p). Following 



Fowler fc Nordheim I (|1928l) . we assume that the triangular shape surface potential is given by V(z) = C — Fz, which is 



changing linearly with ^-coordinate. Here C is the constant surface barrier and F (absorbing the magnitude of electron charge 
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e, we replace eF by F) is the driving field for electron emission from the poles. Consideration of linear type surface barrier 
potential has no-doubt some historical importance. This type of potential barrier was first used in the original work by Fowler 
and Nordheim. However, there are two other important reasons: since this is the first time the problem is solved in presence 
of strong quantizing magnetic field, therefore to get an analytical solution, we have considered such simplest triangular type 
potential barrier. Our intension is to solve the problem analytically in a way in which the physical meaning of the problem is 
not lost. The other reason behind such a choice is that, in the case of strongly magnetized rotating neutron stars / magnetars, 
the produced electric field at the polar region is approximately constant for uniform magnetic field strength at that region 
and constant rotational period of the object. We believe that for the potential barrier at the poles both assumptions are 
approximately valid. Then assuming the conservative force field relation dV/dz = — eF we get the triangular type potential 
barrier at the poles. Under such situation, the Schrodinger equation satisfied by the electrons which are confined within the 
matter (in this case within the neutron star/magnetar crustal matter) is given by (throughout the paper for the sake of 
convenience we assume natural units, i.e., % = c = 1) 
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(1) 



Whereas for the electrons just liberated out through quantum mechanical tunneling, one has to consider the potential V(z) 
along with E. The relevant equation is given by 
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-E + V(z) V = 



(la) 



where the energy eigen value E is given by the eqns.(4) and (5) for two different physical situations. If we assume a separable 
solution for the eqns.(l) and (la), satisfied by freely moving electrons and moving under the potential V(z), the wave functions 
can be represented by 



ip(p,9,z) = 4> n „,m{p,0)f u (z) 



(2) 



where for eqn.(l), the longitudinal part is plane wave type, whereas for eqn.(la), we shall evaluate f v {z) using the technique 
as discussed below. Since there is no potential associated with the transverse motion for the electrons, the free transverse part 
of the wave function is given by 



Me) = 



exp(im#) _!_ 

(27T)V2 P ° 



m +n„ 



2l m U 
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exp 
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2Pc 2 , 



(3) 



with po = (2/eB) 1 / 2 is the Larmor radius, the radius of the lowest Landau orbit and L n m is the Laguerre polynomial. The 
energy eigen value is given by 

„2 



E = -p- + p B B{2n p + v ± v + 1) 
2m e 

without the electron spin and with the inclusion of the electron spin, it will be 



(4) 



E 



2m, 



+ n B B(2n p +u±v + l)T PbB 



(5) 



where ps = e/2m, the Bohr magneton. For the motion of electrons along z-direction, the Schrodinger equation satisfied by 
fv(z) with and without potential V(z) can be obtained by averaging eqns.(l) and (la) respectively over the transverse wave 
function (f>n p ,m(p,d)- Then for just tunneled out electrons in the lowest Landau level with n p = v = m = 0, moving along 
z-axis under the influence of the potential V(z) and having no spin contribution, we have 

1 d 2 /o 



+ V(z)f = {E- pbB) f = wf 



(6) 



2m dz 2 

Whereas for the same kind of liberated electrons with spin polarization in the negative direction of z-axis, we have w — E 
and for the polarization along the direction of magnetic field, w — E — 2psB. On the other hand, for the electrons confined 
within the crustal matter of magnetars and moving freely, we have V(z)=0. The corresponding Schrodinger equation is given 
by 

^+^/o = (7) 

Where Wk = (2mm) 1 ' 2 , is the equi valent electron mome n tum a long z-axis and w can have three possible values as mentioned 
above. Then following eqn.(10) of iFowler fc Nordheim I (|1928l), we can write down the solution of eqn.(7) for free electrons 
inside the crustal matter in the form 
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fo 



1/2 



[oexp(iwfcz) + a exp(— iw^z) 



(8) 



where a is the probability amplitude for electrons moving along the positive direction of z-axis (incident part), whereas a' is 
the corresponding quantity for left moving waves (reflected part from the surface barrier). Next for the electrons just tunneled 
out, we make the following coordinate transformation 

C 



V = 



/ C - w \ 
{-— + Z ) 



(2m e F) 



1/3 



Then we have from eqn.(6) 

d 2 fo 
dy 



2 +y/o = o 



(9) 



(10) 



Now following eqn.(8) of Fowler fc Nordheim ( 19281 ) the solution for this equation is given by 

fo(y) = y 1/2 H[%(ly^) (11) 

where H, 2 ,{(x) is the Hankel function of second kind of order 1/3 with argument x. As we shall see below the quantity Q 
defined in lFowler fc Nordheim I (jl928T ) after eqn.(12) is much larger in our formalism in presence of strong quantizing magnetic 
field, unless the direction of spin polarization for the emitted electrons are opposite to the direction of external magnetic field. 
In absence of electron spin term or polarization along the direction of magnetic field, the factor Q is virtually infinitely large 
(Q « oo). We shall further show that this infinitely large Q will make the transmission probability of electrons vanishingly 
small. In our formalism the factor Q is defined as 



= | (2m, 



F) 1/2 ( 



C 



F 



3/2 



and is related to the Hankel function argument by the relation 



H 



(2) 
1/3 



exp 



at z 



(11a) 



(116). 



Now the transmission probability for electrons as defined in Fowler fc Nordheim (1928) is given by 



D(w) = 

Since the transmission coefficient D related to Q factor by the relation D 
electrons from the zeroth Landau level is related to D by the integral 

eB 
2^2 



(12) 

exp(— 2Q), and the field emission current for the 



(12a) 



with f(w) the Fermi distribution function, then it is quite obvious that in the electron field emission current, the polarization 
effect will come only from the transmission coefficient D(w). In the following, we therefore focus our study in the investigation 
of the properties of the emission coefficient. For large Q, we have from eqn.(18) of Fowler fc Nordheim ( 19281 ) in our modified 
form 



D(w) 



[w{C-w)] 1/2 



C 



■ exp 



. {2me F) [—^-) j 



[w(C- 



^/2l 



c 



■exp(-2Q) 



(13) 



Let us now analyze the argument part of the exponential. Following IFowler fc Nordheim 

I dl928h. we put C = He + Wf, w here 



/i e is the electron Fermi energy and Wf = w c x (B/Bc) 1 ^ in eV is the work function (see lGhosh fc Chakrabarty (2011)) for 



] Jl928; 
;e lGhog 

W CV! A 



the emission of electrons along the direction of magnetic field, where w c w 82.93 and Be w 4.43 x 10 G, the typical value 
of magnetic field strength at which the Landau levels for the electrons are populated in the relativistic scenario. Again, as 
defined before, the quantity w can have three possible values. To get an order of magnitude estimate for the terms containing 
Bohr magneton and work function, we assume that the emitted electrons carry the maximum possible energy, i.e., the electron 
Fermi energy for temperature T — > limit. Then C — w = Wf+l3x {ibB, where the parameter ft = 0, or = 1 or = 2 for 
the spin polarization opposite to the direction of external magnetic field, i.e., with conventional direction of spin polarization 
in presence of strong magnetic field, or no spin term or spin polarization along the direction of magnetic field respectively. 
Taking into account the denominator F of the argument and expressing in terms of magnetic field strength as defined at the 
beginning of the introduction and using /is ~ 5.79 x 10 _15 MeV G _1 , the numerical value for the Bohr magneton, we have 
approximately from the expression for Q as defined above 



h 1 / 2 ^ x 0.5 + w c h] 



1/2 



x nr 6 ) 3/2 x W = q b + q 



(14) 
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where hf = B/B^\ Qb and Q w f are the contributions from Bohr magneton or spin term and work function part respectively. 
It is quite obvious that the contribution from spin term Qb is extremely large for ft — 1 or j3 — 2 and makes the transmission 
coefficient exactly zero. Which physically means that if we do not consider electron spin or assume electron spin polarization 
along the direction of magnetic field, the electron transmission coefficient and in turn the electron transmission current vanishes 
exactly, i.e., there will be no field emission under such situations. On the other hand, the second term, the work function 
part, unlike the first term, gives finite contribution to transmission coefficient. The first term, i.e., the spin term will make Q 
extremely high, even if we do not assume that the emitted electrons are at the top of their Fermi level. This particular term 
is therefore making the cold emission probability of electrons from the poles of strongly magnetized neutron stars exactly 
zero for the non-zero values of /?. For the transmission of spin polarized electrons j3 — and consequently Qb = 0. Therefore 
Q = Qwf - In this situation the electron energy eigen value obtained from eqn.(5) for n p = v — m — is given by 

E = -^- (15) 
2m e v ; 

In this case there will be electron field emission with their spins polarized opposite to the direction of external magnetic 
field. From the above expressions for electron energy eigen value, it is quite obvious that for our model on cold emission of 
polarized electrons the rest of the mathematical formulation will be almost identical with that of Fowler and Nordheim in 
presence of strong quantizing magnetic field. We have further noticed that the values for transmission coefficient D and the 
transmission current R remain finite and large enough for the magnetic field strength < 10 15 G. Of course in our model the 
magnetic field dependency of Q will come from the work function w / and the field intensity F and only the polarized electrons 
are allowed to tunnel through the surface barrier. Then following Fowler and Nordheim we have obtained the cold electron 
emission current from eqn.(12a) at T — > using Sommerfeld' s lemm a. While calculating emission current numerically, we 
have used the exact form of D as given in iFowler fc Nordheim I diiiih after eqn.(16). In fig.(l) we have shown schematically 
the variation of cold electron current with the strength of electric field intensity at the poles and is represented by the solid 
curve. Since the electric field at the poles is produced by the rotating magnetic field, we have also shown in the same graph 
with dashed curve the variation of electron cold current with the intensity of magnetic field. Since the electric field intensity 
varies linearly with the magnetic field strength, the qualitative nature of the curves are almost identical. For the neutron stars 
with very low surface magnetic field, the produced electric field, which acts as driving force, is also small enough, therefore 
the electron field emission current will be extremely small as shown in the curve. For the neutron stars with moderate surface 
magnetic field strength (B ~ 10 12 — 10 15 G) the field current is quite high. Now for the objects with large surface magnetic 
field (B 3> 10 15 G), the work function will also become large. Therefor beyond some maximum value for electron field current, 
since the work function part dominates over the driving electric force at the poles, the electron field current will decrease 
and will become vanishingly small. This is the case for the objects with ultra- high surface magnetic field (3> 10 15 G). In the 
figure, we have not shown the variation of electron field current with the magnetic field strength (computed at the surface) 
for a particular neutron star. Therefore each magnetic field / electric field points corresponds a particular type of compact 
magnetized object, with surface magnetic field from very low to ultra-high values. From the curves it is quite obvious that for 
magnetars with surface magnetic field ~ 10 15 G, the electron field current is quite high and very close to the peak value. Unlike 
the original work of Fowler and Nordheim the tunneling coefficient does not follow exponential law. However, the variation of 
cold current with the electric field strength can be obtained numerically. The numerically fitted functional form is given by 

R = O.26.F3/ 2 exp(-9.8 x F14) (16) 

where R is the field current in Amp/cm 2 , Fu = 10~ 14 F and F 24 = 10~ 24 F. 



3 EFFECT OF STRONG QUANTIZING MAGNETIC FIELD ON COLD EMISSION: RELATIVISTIC 
SCENARIO 

In the relativistic scenario we have repeated the non-relativistic calculation for the cold emission transmission co-efficient. In 
this section we have considered emission of high energy electrons from the polar region of strongly magnetized neutron stars 
with magnetic field 10 12 < B < 10 17 in Gauss. The potential is introduced by hand in the Dirac equation using standard 
relativistic hadro-dynamic technique. Following Fowler and Nordheim, here also we have considered a triangular type potential 
barrier at the polar region. Like the previous case, we have considered cylindrical coordinate system and the choice of gauge 
is same for A. The radial part of upper component satisfies the equation 



d 2 19 ,22 
ap z p op 



R{p) = (17) 



where fi\ — E 2 — m* 2 — p 2 + 2Xk with A = ±1 for up and down spin states respectively, k = eB/2, E is the energy eigen 
value for the Dirac equation and m* — m + V(z) — m + C — Fz is the effective electron mass and m* = m for free electrons 
when V(z) — 0. The solution of the above equation is given by 
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R(p) = Nexp(-~^L n (t) (18) 

where L„(t) is the Legendre polynomial of order n, t = k 2 p 2 and TV is the normalization constant. For up and down spin 
states, the energy eigen values are given by Ef = [p 2 + m 2 + 2neB] 1 ^ 2 and E± = [p 2 + m 2 + 2{n + l)e_B] 1//2 respectively. In 
general we may write E„ = (p 2 + m 2 + 21/eB) 1 ^ 2 . Then following the same averaging technique as discussed in the previous 
section for the non-relativistic case, the wave function of electrons corresponding to the motion along z-direction is given by 



d 2 f„ 



+ (E 2 - m* 2 - 2veB)f v = (19) 



dz 2 

Using the transformation 

m + C -uF l/2 
Z= F 

with the new variable being u, the above equation for f u reduces to 



(20) 



d 2 j v 



+ (a 2 -u)f v =0 (21) 



du 2 

where a = (E 2 — 2feB) 1 ^ 2 / F 1 ^ 2 . which is the well known form of differential equation for one dimensional quantum mechanical 
harmonic oscillator. With a 2 = 21, we have E 2 — 2(yeB + IF), and the solution is given by 

■i 2 
2 



/„,, = iiVexp ( -\ J Hiiu) = (J v )j (say) (22) 



where iV is the normalization constant and Hi{u) is the Hermite polynomial of order I. This spinor solution is for those 
electrons which have already been liberated out through the surface into vacuum under the influence of electric field F (here 
liberated out from the crustal matter of strongly magnetized neutron stars or magnetars into the magnetosphere through 
polar region). 

The equation satisfied by free electrons bound within the crustal matter by the barrier potential at the surface can be 
obtained by putting V(z) — and is given by 

where in t his case, a = (E 2 — m 2 — 2 1/eB) 1 ^ 2 is the free electron momentum along z-axis of energy E. Now following the 
notation of Fowler &: Nordheim ( 1928h . we express the solution for free electrons within the system, confined by the surface 



barrier V(z), in the form 

f" = ^T72 [ aex P( iaz ) + a'exp(-zaz)] = (/„)„ (say) (24) 

where as before a is the probability amplitude for electrons moving along the positive direction of 2-axis (incident part), whereas 
a' is the corresponding quantity for left moving waves (reflected part from the surface barrier). Assuming the interface between 
the crustal matter of the strongly magn etized neutron stars and th e magnetosphere is at z = 0, the wave function and their 
derivatives must be continuous at z = o lFowler fc Nordheim ( 19281 ). i.e., 



(/*(0))„ = (MO)), and {fl{G)) II = {fu{Q)) I (25) 
Using the relation H[{u) = 2fi?;_i(u), we have 

a + a' = Na 1/2 expf-^\H,(u ) (26) 
and 

ia 1/2 (a - a') = iVexp (~ [«ofli(uo) - 2lHi- 1 {u )\ (27) 
where Uo = u(z = 0) = m* /F 1 ^ 2 . These two conditions may be rearranged in the following form 

a + a' = X and a — a — iY (28) 

where X and Y are the two real quantities. Hence it is straight forward to verify from eqn.(12) that in the relativistic scenario 
the transmission coefficient vanishes exactly. Therefore from the analysis of this section, we may conclude that if the barrier 
in combination with the external electrostatic driving force behaves like a scalar type potential and is triangular in shape at 
the surface, then the relativistic electrons can not tunnel through the surface barrier whatever be their kinetic energies and 
the strength of external electric field. 
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4 CONCLUSION 

The non-relativistic scenario of cold electron emission in presence of strong magnetic field is believed to be the first attempt 
in this direction. While obtaining the electron transmission probability in the non-relativistic scenario under the influence of 
strong electric field at the poles, we have noticed that in our theoretical formalism, the emission is allowed if we take electron 
spin into account and also the electrons have conventional spin polarization, i.e., opposite to the direction of external magnetic 
field. Even if all these criteria are satisfied, at extremely high magnetic field strength since the electron work function becomes 
large enough, the transmission coefficient drops to zero. Since there are no such stellar objects with surface magnetic field 
strength S> 10 5 G, the vanishingly small transmission coefficients will therefore not be possible in reality. The low charge 
density magneto-sphere will therefore only be possible if such super exotic compact stellar objects with ultra-strong surface 
magnetic field exist in nature. In the conventional magnetar or strongly magnetized neutron star case, the electron field 
current is quite high, very close to the peak value, which is also obvious from the figure. As a result there will be enough 
curvature gamma photons produced by the energetic electrons, which in turn produce enormous amount of e~ — e + pairs 
to form normal magneto-spheric plasma. In our model the only difference from non-magnetic or low magnetic emission case 
is that the primary electrons are spin polarized in the direction opposite to the direction of external magnetic field. As a 



consequence, for the conventional magnetar case, the results obtain in iBeloborodov et al. I (|2007f ) on the formation of corona 



of magnetars will not be seriously affected. Therefore, in the magnetar magnetosphere of a neutron star with B < 10 5 G, 
the primary electrons are all spin polarized along — B. Since the electron emission current is almost zero for the objects with 
ultra-high magnetic field strength (S> 10 G), then if the electron emission process in such exotic objects is dominated by 
cold emission, the charge density of e~ — e + -plasma in the magnetosphere will be extremely low. As a consequence there 
will be very weak synchrotron emission in the radio wave band. The other possible mechanism by which e~ — e + -plasma can 
be produced in the magnetosphere of magnetars are (i) thermoelectric emission of electrons from the polar region and (ii) 
photo emission from the same region. However the work function at the polar region of a typical magnetar is several GeV, 
whereas the temperature can be at most a few hundred MeV for an young magnetar. The thermionic emission will therefore 
be suppressed by the Boltzmann factor exp(— Wf/kT). Whereas in the case of photo emission, the energy of the induced 
photon must be of GeV order (7-photons). At the polar region, if any such photons exist to liberate electrons, they must have 
produced as curvature photons by high energy electrons moving along the magnetic lines of forces. Since the possibility of 
such electrons is very rare, the number of high frequency photon is also vanishingly small, as a result there will be almost no 
creation of secondary e~ — e + -pairs in the photo emission process. 

In conclusions, from our relativistic formalism of cold emission of electrons we can state that relativistic electrons popu- 
lating the neutron star interior can not be extracted from cold emission from the poles of a neutron star, independently from 
the magnetic field strength. Non-relativistic electrons with anti-parallel spin can be extracted for standard (observed) values 
of magnetic field strengths, but can not be extracted from the surface of objects with B ^> 10 5 G (in the eventuality that 
such exotic objects can exist). 
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Figure 1. (i) With dashed curve the variation of electron field current plotted along left side vertical axis with the neutron star magnetic 
field plotted along x-axis at the bottom and expressed in terms of critical field for electron is shown, (ii) The variation of same 

quantity plotted along right side vertical axis with the electric field in Volt/cm produced by rotating magnetic field of magnetars, plotted 
along horizontal axis at the top is shown by solid curve. 
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